High capacity spatial multimode quantum memories based on atomic ensembles 
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We study spatial multimode quantum memories based on light storage in extended ensembles of 
A-type atoms. We show that such quantum light-matter interfaces allow for highly efficient storage 
of many spatial modes. In particular, forward operating memories possess excellent scaling with 
the important physical parameters: quadratic scaling with the Fresnel number and even cubic with 
the optical depth of the atomic ensemble. Thus, the simultaneous use of both the longitudinal and 
transverse shape of the stored spin wave modes constitutes a valuable and so far overlooked resource 
for multimode quantum memories. 
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Introduction: Photons are ideal candidates for carry- 
ing quantum information, however, in order to process 
the information or store it for a later time, a quantum 
storage medium is needed. To achieve this, one needs to 
establish a controllable and efficient light-matter inter- 
face that will store light as a stationary excitation in the 
atomic medium while preserving quantum correlations. 
Quantum memories have already been demonstrated in 
a number of experiments based on atomic ensembles, see 
e.g. [ll-Q as well as solid state systems, see e.g. 0- 
0]. Most of the realized memories support only a sin- 
gle mode, but in order to constitute a building block 
for a "quantum internet" and for quantum compu- 
tation they should preferably be capable of storing as 
many modes as possible ll|, [l^ . The protocols for mul- 
timode memories can be based on various degrees of free- 
dom: spatial or directional [31 modes as well as on 
controlled reversible inhomogeneous broadening (CRIB) 
isl Til and on atomic frequency combs (AFCs) flTll . the 
latter of which has been realized experimentally |l8| . 

In practice, experimentalists have access to limited 
physical resources and typically work with a given atomic 
ensemble where one wants to make full use of the capa- 
bilities that the available physical system has to offer. 
Limitations to the efficient storage of many modes can 
result, for instance, from the achievable density of cold 
atoms, imperfection of the optical pumping or the geom- 
etry of the atomic cloud. Typically in the experiments 
a single or a few independently acting transverse modes 
are employed. 

Here, we study how to make full use of the additional 
resource given by the spatial extent of the atomic ensem- 
bles. We show that the simultaneous use of the longitu- 
dinal and transverse degrees of freedom allows for highly 
efficient storage of many spatial light modes resulting in 
capacities higher than previously expected. The number 
of modes one can store with high efficiency depends on 
the choice of the direction of retrieval relative to that of 
the storage process. We demonstrate that the forward 
operating memory, with the retrieved light traveling in 
the direction of the input signal, provides an excellent 



multimode memory resource. It has a remarkable scal- 
ing with the important physical parameters: the optical 
depth do and the Fresnel number of the atomic ensem- 
ble F. Earlier estimations of the capacity of multimode 
memories predict a scaling with ^fdo for Raman mem- 
ories as well as for protocols based on electromagneti- 
cally induced transparency (KIT) and a linear scaling for 
the CRIB protocols [l^. The dependence on the Fresnel 
number has only been roughly estimated in Refs. (l2l.[20j 
to be ~ . Here, we show by direct calculation that 
the simultaneous use of the transverse and longitudinal 
shape of the stored spin wave mode leads to quantum 
memories with capacities that scale quadratically with 
the Fresnel number (~ and even cubically with the 
optical depth of the atomic ensemble (~ dp). This is 
much better than expected from previous results. For 
comparison, we also study the backward operating spa- 
tial memory, with the retrieved light traveling in the op- 
posite direction of the input light, and show that it can 
also serve as a high capacity multimode memory but with 
a slightly less promising scaling with the physical param- 
eters. 

Model: In order to analyze the capacity of the spa- 
tial multimode quantum memory, we use the three- 
dimensional theory for quantum memories based on A- 
type atomic ensembles presented in Ref. 2l|. There, it 
was shown that the crucial physical parameters deter- 
mining the quality of the quantum memory are the op- 
tical depth do and the Fresnel number of the atomic en- 
semble F. We consider a cylindrically symmetric atomic 
ensemble with a Gaussian distribution in the radial direc- 
tion n{p) = riQ exp[— p^/(2fT^)] [see Fig. [1] (left)], where 
no — Np^/{2-KLa\), Na is the number of atoms, L is the 
length and a-± describes the width of the ensemble. The 
density along the longitudinal z axis has been assumed 
constant for simplicity. The geometry of the ensemble is 
described by its Fresnel number F = ^'^/{XqL), where Ao 
is the wavelength of the light. The weak quantum field 
carries the quantum information to be stored into the 
atomic ensemble and couples states |0) and |e) with cou- 
pling strength 5 [see Fig. [1] (right)]. States |1) and |e) are 
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In order to solve the equations of motion ([Ij, we 
perform the Laplace transform in time £{g{t)} = 



roc 

Jo 



*g{t)dt^ which allows us to eliminate the differen- 



tial equation for light. Here, we have assumed a constant 
driving field in space and time f2. In consequence, we 
can write the relation between the light modes and the 
spin wave (here for the forward operating memory) in the 
form of input-output beam splitter relations 



FIG. 1. (Color online) (left) Atomic ensemble of cylindri- 
cal symmetry with Gaussian atom distribution, (right) A 
schematic plot of a level scheme of the considered A-type 
atoms. 



coupled by the strong classical control field, which sets 
the propagation direction of the retrieved light. VL{t) is 
the Rabi frequency of the driving field and A denotes the 
detuning from the excited state |e), which spontaneously 
decays at a rate 7. The information is stored in a collec- 
tive state of the ensemble represented by stationary spin 
wave excitations described by 5 ~ |0)i(l|. 

For simplicity, we solve the three-dimensional problem 
of the multimode quantum light-matter interface within 
the adiabatic approximation, where the excited state |e) 
is eliminated. The equations of motion within the parax- 
ial approximation (in the co-moving frame t' = t — z /c) 
read [2l| 
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Here, we have introduced the dimensionless time t = "ft' , 
detuning A = A/7, position z = z/L and Rabi frequency 
f2(t) — Q{i)/j. The peak optical depth do = 4Lno|gp/7 
quantifies the absorption of the resonant light in the ab- 
sence of the control field U{t). We omit here the quan- 
tum noise since it is not needed for calculating the ef- 
ficiency below ,22]. The slowly varying light d{z,t) and 
spin wave S{z, t) operators have been expanded in a basis 
set of transverse mode functions Umn{ri_) (see Ref. [2lj 
for detailed derivations). Due to the sample symme- 
try as well as for numerical reasons, we chose a set of 
Bessel beams indexed by n and the azimuthal quan- 
tum number m. The last terms in the above equa- 
tions of motion describe the coherent interaction be- 
tween light and matter, which is quantified by the op- 
tical depth, Rabi frequency, detuning and the coupling 
matrix Bmn^m'n' = / d?rx_u*^^{ri^)um,n'{ri^)n{ri_)/nQ. 
Note, that we use here a compact vector and matrix no- 
tation where amn{z^ t) — > a(z,i) and Bmn.m'n' — ^ B- 



flout (t^) = / dzK[Q,,uj,z] So{z), 
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So{S) = TT- / dCoK^p*,^, 1 - i]ai„(d;). (3) 
27ri 

The transformation matrix IK depends on frequency cj, 
position z and the physical parameters of the system 
fl,A,F and do- The analytical expressions for the ma- 
trix as well as beam splitter relations for the backward 
read-out are presented in detail in Ref. [2l|. The effi- 
ciency of the quantum memory is the ratio between the 
number of outgoing and incoming light field excitations, 
which, assuming a normalized incoming light mode, can 
be written as 



df|aout(?)|' 
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Here, the kernel matrix M gives full information about 
the relation between the input and output modes for 
a given set of physical parameters. After discretizing 
frequency and position, we diagonalize the large ker- 
nel matrix, which for the forward retrieval is M[i/, v'] 



z']. A 



set of characteristic efficiencies is thereby obtained as the 
eigenvalues of this matrix together with the correspond- 
ing set of incoming light modes d[n{u). The optimal in- 
coming light mode to store into the atomic ensemble cor- 
responds to the eigenvector with the highest eigenvalue, 
which is the maximal efficiency of the memory. The 
remaining eigenvectors correspond to orthogonal modes 
that can be stored with lower efficiencies. We show in the 
following that, in general, there exist many light modes 
that give high efficiencies of the quantum memory. 

The figure of merit for multimode quantum memories 
is the capacity, which can be defined in at least two dif- 
ferent ways. Firstly, one can simply count the number 
of modes with an efficiency above a minimal value ?7min- 
Secondly, a more sophisticated measure can be obtained 
from the quantum capacity of a Gaussian channel with ef- 
ficiency 77, Q(ry) = max{0,log2 |?7|-log2 |1-??|} Qiv) 
is the average number of qubits that can be perfectly 
stored and retrieved from a particular mode with com- 
bined storage and retrieval efficiency 77, provided many 
copies of such memory and optimal encoding, decoding 
and error correction of the stored information. The ca- 
pacity of the memory is obtained by summing the capac- 
ities for all modes; only modes stored and retrieved with 
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Fresnel number, F 



FIG. 2. (Color online) The number of modes with efficiency 
77 > 0.5 (diamonds and stars) and the capacity Cm (squares 
and circles) for different values of the azimuthal quantum 
number m for forward (blue, bottom and left axis) and back- 
ward (black, top and right axis) operating memory with 
do = 100 and F = 1. 



a combined efficiency above r] = 0.5 contribute due to 
the no-cloning theorem [24], so that C — X]i;>o bQiv)- 

Results: In order to gain insight into the multiniode 
character of the spatial quantum memory under consid- 
eration, we first calculate the capacity and the number 
of good orthogonal modes within subspaces with a fixed 
azimuthal quantum number m, see Fig. [2j All numerical 
results are obtained for resonant memories with A = 0. 
Note, that we show only the non-negative values of the 
azimuthal quantum number m since the modes for m 
and — m have the same efficiencies. The number of modes 
with high efficiency decreases with growing \m\ since light 
beams with |m| > vanish to increasing degree towards 
the center of the atomic ensemble. Therefore, it is harder 
to focus the input light mode into the dense center of 
the atomic cloud, leading to a decreased effective optical 
depth do. One can see that for do = 100 and F = 1, the 
forward operating memory yields higher capacities than 
the backward one. This can be explained by the fact 
that the driving light used for the backward retrieval re- 
verses the longitudinal phase of the stored excitation but 
cannot properly reverse its transverse profile except for a 
spin wave with a uniform transverse phase. Thus in the 
case of any transverse phase gradient of the stored spin 
wave, the irreversible transverse phase leads to unwanted 
diffraction effects f2l|, which in consequence reduce the 
efficiencies and the number of modes that can be stored 
in the backward operating memory. 

The total capacity of the forward operating quantum 
memory Cf is presented in Fig. [3] (red closed circles 
and stars) as a function of the Fresnel number of the 
atomic ensemble F for two values of the peak optical 
depth, do = 40 and 100. In both cases, the total capac- 
ity Cf reaches high values and grows quadratically with 
the Fresnel number, Ct ~ F^. To investigate whether 



FIG. 3. (Color online) The capacity Cf for forward (red 
closed circles and stars) and Cb backward (black open cir- 
cles and stars) operating memories as functions of the Fresnel 
number F for do ~ 40 (circles) and do = 100 (stars). The 
number of modes for the forward memory A^m.f for do = 100 
with the threshold efficiency ?7min = 0.5 (light gray stars) and 
'?min = 0.6 (light gray triangles). The dashed lines are the 
quadratic Cf ~ and linear Cb ~ -F fits for the forward and 
backward direction of the read-out, respectively. 

the scaling is independent of the two capacity measures 
mentioned here, we also plot the number of modes Nf 
for do = 100 for two values of the threshold efficiency 
'7min = 0.5 and 0.6. We see that the quadratic scal- 
ing with the Fresnel number F is universal and either of 
these N or C may be used as the appropriate measure of 
the capacity. 

The scaling with the Fresnel number depends, how- 
ever, on the direction of the read-out. To compare the 
two operating modes of the memory, we also calculated 
the capacity of the backward memory (black open circles 
and stars). We find that in this case, the capacity only 
scales linearly with the Fresnel number but also achieves 
high values. Thus, even though the highest single mode 
efficiency for larger Fresnel numbers is achieved for back- 
ward operating memories {2l'| , the number of good modes 
is larger for the forward memories in this F regime. 

The scaling of the capacity can be explained by consid- 
ering the diffraction of light in the atomic medium. For 
the forward operating memory, the divergence angle of 
the incoming light beam is 9 ^ Xo/rj_, where r±_ is the 
transverse waist of the stored stationary excitation. The 
mode can in this case be so focused that the maximum 
divergence angle becomes limited by the geometry of the 
ensemble: tan^max ~ ^max ~ '^±/L. From this we obtain 
that the minimal achievable waist of the stored excita- 
tion is r^ inin ^ XoL/a±. In consequence, the capacity of 
the memory, which is proportional to the ratio between 
the cross section area of the ensemble and the minimal 
waist squared, Cf ^ '^'±/^±mm — ^ leading to the 
quadratic dependence on the Fresnel number. As men- 
tioned above, in the case of backward operating memory, 
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FIG. 4. (Color online) The number of modes for the for- 
ward operating memory with efficiency 77 > 0.65 for two 
values of the Fresnel number F = 0.2 (black triangles) and 
F — 0.5 (green stars) and the capacity of the backward op- 
erating memory for F = 0.2 (red circles) and F = 0.5 (blue 
squares). The dashed lines are the corresponding fits for for- 



ward N{ do and backward Ch 
out, respectively. 
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dn direction of the read- 



the problem of irreversible transverse phase arises [21 1. 
Therefore one cannot focus the light beam as strongly 
as in the case of the forward operating memory. Requir- 
ing the phase to be constant across the transverse profile 
leads to 0,nax ^ r±/L and in consequence to the capacity 
of the memory Cb ^ cr^^/^'l = F, showing the linear 
dependence on F. 

Above, we have presented the multimode character of 
the transverse degrees of freedom in the quantum mem- 
ory and its dependence on the Fresnel number of the 
atomic cloud. Now, we show that including simulta- 
neously the transverse and longitudinal modes leads to 
much stronger scaling of the capacity with the optical 
depth do than expected from considering these two de- 
grees of freedom separately. We have optimized the full 
three-dimensional quantum memory and found that the 
capacity of the forward operating memory has a promis- 
ing cubic scaling with the optical depth C{ ^ dp, see 
Fig. |4l much better than previously expected. (The rea- 
son for using the number of modes A^f with an efficiency 
of at least r/min = 0.65 in the figure was purely numeri- 
cal. Large values of F and do make it hard to calculate 
all modes with 77 > 77,„in = 0.5 as required for the total 
capacity C.) 

This remarkable cubic scaling of the capacity of the 
forward operating memory with the optical depth do can 
be understood by means of the phase space volume of 
the incoming light mode Vps ~ A±,max?'max, where k± 
is the maximal perpendicular wave number and r 
the maximal waist of the incoming light beam. If we 
consider light modes coming with different angles with 
respect to the ensemble axis z, the corresponding effec- 
tive optical depth for a given angle can be expressed 
as do{d) ~ docr±/{L0) for 1 > 6* > a±/L. Storage 
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with a certain efficiency requires at least some critical 
optical depth dcrit, which leads to a maximal permit- 
ted divergence angle 6'max ~ d ^do / {Ldcrit) ■ Therefore, 
the maximal perpendicular wave number is proportional 
to the optical depth, fcj.,max ~ crj_(io/(LAodcrit), which 
we have confirmed numerically. For the same reason, 
the maximal waist of the incoming light beam r^ax ^ 
tan^max-^ ~ Oma-xL ^ cr±d/ dcrit also scales linearly with 
the optical depth. In consequence, the number of avail- 
able modes for each azimuthal quantum number m is pro- 
portional to the phase space volume Vps ~ Fd^/d^j.^^. On 
the other hand, the maximal azimuthal quantum num- 
ber TOmax ~ o'_LA;_L_max ~ Fdo/dait scalcs linearly with 
the optical depth as well (also confirmed numerically). 
Finally, we find that the total number of good modes, 
expressing the capacity of the forward operating mem- 
ory under consideration, Nf ~ WmaxV^s ~ ^^'^o/'^CTit 
scales quadratically with the Fresnel number F and cu- 
bically with the optical depth do. This is a much bet- 
ter scaling than the one achievable for other types of 
multimode memories with AFC or CRIB, where a lin- 
ear scaling with the optical depth is obtained [l9|. For 
comparison, we also provide the results for the backward 
operating memory, see Fig. HI Here, the scaling is less 



but still high values of the capacity 



3/2 

promising Cb ~ do ' 
are achievable. 

Conclusion: We have calculated the capacity of spatial 
quantum memories based on A-type atomic ensembles 
and thereby shown that they allow for storage of many 
light modes and exhibit a remarkable scaling with the 
important physical parameters. For memories operated 
in the forward direction, the capacity scales quadratically 
with the Fresnel number F and cubically with the opti- 
cal depth of the atomic ensemble do , which is much bet- 
ter than previously expected These results reveal 
that the transverse degrees of freedom combined with 
the longitudinal ones constitute a valuable resource for 
multimode quantum memories with excellent capacities. 
These results can be directly used in current experiments 
with extended ensembles of A-type atoms. 
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